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Dark Energy is some of the weirdest and most mysterious stuff in tlie universe tliat tends 
to increase the rate of expansion of the universe. Two commonly known forms of dark energy 
are the cosmological constant, a constant energy density filling space homogeneously, and 
scalar fields such as quintessence or moduli whose energy density can vary with time. We 
explore one particular model for dynamic dark energy; quintessence driven by a scalar dilaton 
field. We propose an ansatz for the form of the dilaton field, |(/)(a)|/mp = a^^ Int + a^t" = 
alna + /3a^^, where a is the scale factor and a and are parameters of the model. This 
phenomcnological ansatz for can be motivated by generic solutions of a scalar dilaton 
field in many effective string theory and string-inspired gravity models in four dimensions. 
Using a compilation of current data including type la supernovac, we impose observational 
constraints on the slope parameters like a and C and then discuss the relation of our results to 
analytical constrains on various cosmological parameters, including the dark energy equation 
of state. Sensible constraints are imposed on model parameters like a and C as well as on the 
dark energy /dark matter couplings using results from structure formation. The constraints 
of this model are shown to encompass the cosmological constant limit within Icr error bars. 

PACS numbers: pacs: ll.25.-w, 04.80. Cc, 98.80.Cq 



I. INTRODUCTION 

The discovery of a cosmic acceleration of the universe {l-4] has added a new challenge for 
fundamental theories of physics and cosmology. NASA's observations ^7] show that the kind 
of matter of which stars and galaxies are made forms less than 5% of the universe's total mass. 
Several independent observations indicate that about 73% of the total energy density of the 
universe is in the form of mysterious dark energy or gravitationally repulsive energy, and about 
22% is in the form of non-baryonic cold dark matter particles which clump gravitationally, but 
which have never been directly detected. These scientific enigmas suggest we should look to new 
physics beyond the theories of the standard model of particle physics and Einstein's theory of 
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General Relativity. 

The cosmic acceleration of the universe is often attributed to the existence of some form 
of dark energy. Understanding the nature of this dark energy is one of the most challenging 
theoretical problems facing present-day cosmology. There is no shortage of ideas for what 
dark energy might be, from the quantum vacuum to a (ultra-) light scalar field, e.g. cor ning 



from modified f(R) gravity scalar-tensor theories 

[2^) or space-time fluctuations of a holographic origin 
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141] ■ string-inspired cosmologies 
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- 24( 1 and the likes. In some other 
approaches, cosmic acceleration is a manifestation of new gravitational physics, whi ch may 
involve modification of gravity at long distances 



or warped extra dimensions 



3d, 



25[, cosmology with extra dimensions j26l429l| 



311 ]. There are also proposals which explain the supernova 
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331 ]. The ma 
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ority of 



381 ] ■ has 



observations using gravitational effects of density inhomogeneities 
work in the subject of cosmic acceleration, as well as the reviews of dark energy 
shown that any viable theory of gravity must look in the far infrared, or on astrophysical scales, 
like a combination of scalar fields and gauge fields (weakly) interacting with General Relativity. 
Moreover, since the cosmological backgrounds of interest are time-dependent, for the consistency 
of any dynamical dark energy model or theory of cosmic acceleration, such scalar fields must 



roll but only slowly as time progresses (see Ref. [39| for some related discussion). 

The ACDM model is the simplest known cosmological model for concurrent universe where 
it is assumed that most matter is in the form of cold dark matter, and dark energy is simply 
described by the vacuum energy < Tni, >var^ ^9uu in Einstein's field equations 
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SttG (TJ^atter + ^vac) • 



fl) 



4Ch42|]. so it may be chosen to match 



The value of the constant A is not set by the theory 
with observations. Vacuum energy does not vary with space or time and is not dynamical. 
This solution seems to be the best-fit to most observations, especially in regard the dark energy 
equation of state, but it does lead to a couple of extremely unlikely scenarios leading us to 
question the validity of using a cosmological constant to describe dark energy. The first scenario 
is the "Why now?" problem. At the present epoch the universe is composed of approximately 
27% matter (both baryons and dark matter) and 73% dark energy. If we lived in a slightly 
earlier or later time these values would be different by several orders of magnitude, so how come 
they are comparable at the current epoch? If the dark energy and matter density were similar 
at some earlier epoch, cosmic acceleration would have began earlier and there would not have 
been enough time to form structures, including our galaxy. 

It is intuitive to think of the cosmological constant as a vacuum energy or the energy contained 
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in free space. As the universe expands, more space is created so there is more gravitational 
vacuum energy or free space. As a result the energy density due to a cosmological constant 
remains the same throughout time. But, since a quantum field must be quantized at each 
and every point in space, it is also useful to think of a vacuum as not just empty space but 
containing a simple harmonic oscillator at each point. Any vibrations in this field space can 
propagate through space as a wave. In the absence of ordinary matter, the hypothetical simple 
harmonic oscillators will be in the ground state which has non-zero vacuum energy. Integrating 
the ground state energy over all space and renormalising it gives us the vacuum energy density. 
The theoretical value predicted in this way is, however, larger by a huge factor of at least 
10^^ than the mass scale associated with the observed value of gravitational vacuum energy 
density, which is ~ 10"^ eV. Although the quantity may conceivably be linked with 
small numbers such as the neutrino mass {rriy ~ P^l^)^ or the ratio PEw/'^Planck ~ /'A' ^ 
theoretical prediction for the value of ~ 10"'^ eV is currently lacking. 

The theory confronting a cosmological constant as the source of dark energy is generally 
called "quintessence". In models of quintessence, one assumes that dark energy is dynamical 
in nature instead of being a constant in time and space by utilizing a standard scalar field. 
In all modern unified theories including Kaluza-Klein theory, supergravity and superstring, a 
fundamental scalar field called the dilaton exists along with a spin-2 graviton - the messenger of 
gravitational force. In extra-dimensional theories of gravity, such a scalar field can be the radion 
field that encodes the distance between the two branes (or brane-antibrane pairs) in the fifth 
dimension. There can be many such scalar fields in higher dimensional braneworld theories that 
may cause cosmic acceleration if they were to self interact and generate a positive potential. 
A canonical scalar field (j) having large self interacting potential but almost negligibly small 
kinetic term behaves almost as Einstein's cosmological constant, at least, from a viewpoint of 
its equation of state. 

In addition to exhibiting an effect on the rate of expansion of the universe, a scalar field 
quintessence, if it interacts with dark matter, may lead to an observable effect on how the 
large scale structure is growing. Recent observations show that clusters of galaxies may have 
formed as a result of small density fluctuations in the primordial plasma when tiny quantum 
fluctuations were amplified by inflation to form pockets of space that were more dense than the 
average energy density. Once matter begins to dominate, gravitational collapse would amplify 
the tiny inhomogeneities left by inflation and visible structures begin to grow into the already 
formed dark matter halos. In a scenario where dark energy or a scalar quintessence held interacts 
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with dark matter, the rate of growth of the dark matter halos would be affected, and thus the 
formation of visible structures. 

Another important motivation for considering a dynamical scalar field is an essence of one 
or more dynamical energy components in the early universe: so called "inflaton". This inflaton 
field could easily release its potential energy to matter and radiation as the field rolls down to 
its lowest (or metastable) energy state. It is quite plausible that the late time quintessential 



scalar field (together with dilatonic dark matter) is a remnant of the inflaton field (see, e.g. 431 ]) 
although there is no guarantee they are the same. 

In this paper we demand a specific form for the scalar field itself, which is indeed an approx- 
imation to some generic solutions for a dilaton in some effective string theory models in four 
dimensions. Our choice is motivated from the fact that a similar ansatz would be useful when 
one is considering the early inflation of the universe. Starting with the standard form of scalar- 
tensor theory of gravity, we will derive fleld equations that depend on both the scalar field and 
the coupling of that field to dark matter. We impose constraints on the parameters by fitting 



the model to observational data from the the cosmic microwave background (CMB) shift 

4^]. Observa- 



baryon acoustic osciUations (BAO) [45|], SNIa (Gold-FlIST sample) [46| and SNLS 
tional results from CMB can be used to tightly constrain many useful cosmological parameters, 
including the parameters of inflation. 

The growth of structure is also expected to depend on the way that dark energy is cou- 
pled to dark matter, granting yet another constraint. Our focus in this work is to explore 
ways to constrain the parameters of the solutions, which include curve flttings to the 
WMAP+BAO+SNIa+SNLS data sets. Also, we will consider how the links between quintessen- 
tial dark energy and inflation may constrain the parameters of the model. Finally we discuss on 
how the results on structure formations constrain the model parameters and the matter fluctu- 
ation growth rate. Using those constraints we show that the solutions of the model encompass 
the cosmological constant limit within la error bars. 

II. QUINTESSENCE: THEORETICAL MOTIVATION 

Einstein formulated gravity as a tensor theory, where the field equations of general relativity 
can be obtained by varying the action 

S = lj yf^d^xR + j K^Cm[^,n,gt,^]V^d^x, (2) 
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where K=mn 1 = VSvrG is the inverse of the Planck mass and G is Newton's gravitational 
constant, '^m denotes the (ordinary plus) matter. The first integral is the gravitational part of 
the action and the second integral is the matter part, and Cm is the matter Lagrangian. Varying 
this action with the metric g^y gives Einstein's field equations 

^{R^.u-^Rg^.u^=-^J^ = T^u, (3) 

where T^jy is the energy momentum tensor describing the density and flux of energy. 

Dirac was the first to introduce a scalar field into Einstein field equations, suggesting that 
the Newton's constant could be viewed as a time-dependent parameter. Thereafter Jordan 
and, independently. Brans and Dicke introduced a monomial scalar field known as Brans-Dicke 



dilaton 



48l | , which also assumes a time varying gravitational constant and incorporates a scalar 



field into Einstein's tensor theory. It was motivated by Mach's principle, that "mass there 
influences inertia here". This was interpreted by Dicke to mean the gravitational constant 
should be a function of the mass distribution of the universe, implying that ^ oc or that a 
scalar field somewhat plays the role of the inverse of Newton's constant, ^ cc ip. Unfortunately, 
the Brans-Dicke dilaton, with the action 



(4) 



which is without a field potential for is unacceptable since a massless scalar field could easily 
create a long range fifth force - which, however, does not seem to exist in nature. One can 
therefore modify the Brans-Dicke type gravitational action as 

S = r^ j V^d^xif{ip)R + C^) + J ^d^x/:„[5^„^',„], (5) 

where £^ is a scalar field Lagrangian in Jordan frame. In the particular picture of Brane-Dicke 
model, one defines /(v?) tx ip. This is a special case of a general scalar-tensor theory in the 
Jordan-Fierz frame, where (p is directly coupled to the Ricci scalar, see, e.g. Refs. 49l. ISOl] . 
By using a canonical transformation, we can write the action (j4|) in the Einstein frame 

5" = ^ (^-^ + ^/^d^x + j Cm[g^lu,'^m]^/-^d'^x. (6) 

The scalar field Lagrangian is assumed to have the standard form 

JC-^ = -lg'"'d^^d,^-V{cP). (7) 

The Einstein frame metric g^i, is now related to the Jordan frame metric g^^ via = A(0)^(7^i,, 
where A{(j)) is a function of (j)- One has 7 = ±1 depending on whether the scalar field is the 
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real or imaginary part of some complex modulus or axio-dilaton field present in some higher 
dimensional theories of gravity. Models similar to the one here were studied before, see e.g. 
Refs. 



51 



52l |. but here we do not constrain the form of the scalar field potential rather the 
evolution of the quintessence scalar field (j). Moreover, we would assume that (/> is a canonical 
scalar field (7 = 1) in Einstein frame and it is coupled to both the ordinary and dark matter 
particles but with different gravitational couplings. 

To generalize the action ^ further, we assume that there is a number of different species 
of matter, including cold dark matter and baryons - say j components - each coupling to 
the metric differently with coupling function, Aj[(j)). As the simplest possibility, we choose 
exponential couplings 

Aj{4))(xe''^^^, (8) 

where, in general, Qj can be a function of (j). The matter Lagrangian is now the sum of all 
components, accounting for all the different sectors: 



3 

We are considering here a scenario in which a purely dark sector interaction exists, resulting 
from a nonminimal coupling of dark matter to a fundamental scalar field or quintessence. Such 
couplings give rise to additional forces on dark matter particles in addition to gravity. 

In the early 1990s, Damour, Gibbons and Gundlach 53| showed that a cosmological model 



with two sectors - cold dark matter with coupling Qc{(t>)^ and normal (baryonic) matter with 
coupling <5fe(0) - satisfies both the weak equivalence principle (WEP) and constraints from Solar 
System observations, provided that 



< 10-^ (10) 

where (f)^ is the present-day cosmological background value of (j). On the other hand, since the 
nature of both dark energy and dark matter are still unknown, there is no physical argument 
that excludes a possible interaction between them. Furthermore, the magnitude of the baryonic 
coupling is constrained from radar time-delay measurements, \Qb\ < 0.032 which limits the cold 
dark matter coupling to IQqdmI < 1- 

In this paper, we specialize to models where the couplings are nonnegligible only in the dark 
sectors, in which case the theory satisfies the WEP. As long as the scalar field is coupled non- 
minimally only to dark matter, the model automatically satisfies the solar system constraint 



56( 1 ■ see also Refs. 



571, 



for some related discussions. 
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Let us assume that the spatial curvature of our universe was erased with inflation in the early 
universe and choose a spatially flat Friedmann-Robertson- Walker metric, 

ds"^ = -dt"^ + a{tfdx^, (11) 

where a{t) is the scale factor, normalized such that a = ao = 1 at present, i.e. at t = tQ. This 
metric assumes the universe is homogeneous and isotropic on largest scales, which is a good 
approximation at large scales and is consistent with CMB observations Q]. For simplicity, we 
also assume that all forms of matter, including dark energy, behave as a perfect fluid. Each form 
of matter may be characterized by the fluid equation of state, Wj = Pi/pi- In particular, the 
cosmological constant may be thought of as a fluid with the equation of state wa = Pa/ Pa 
Other forms of matter or energy are assumed to be perfect fluids as well, for radiation or 
relativistic particles, w^ = — 1/3 and for baryons or cold dark matter w^ — 0. The equation of 
state for a scalar quintessence is then defined as 



which varies with time. 

The variation of the action ([6]) with respect to g^^ and </> gives 



3 r72 , 7 ;2 



At2 



+ 1 ^2 ^ y (0) + (0) ^ = 0, (13) 



2ij + 3i/2^ +2^2_^(^)^^4(^) ^(w,p,) = 0, (14) 



2 



where Hit) = d/a is the Hubble parameter. For the derivation of these equations, see 55l. l56l. 1591] . 
In ([15]), it can be seen that (p couples to the trace of the stress-energy tensor, 

T^. = 5^%. = E(l-3wi)p^. (16) 

i 

For radiation (w^ = 1/3), the stress-energy tensor is traceless, so there is no interaction between 
dark energy and radiation. Equations (|13p -()15p can be supplemented by a fourth equation 
arising from the conservation of energy- momentum tensor, which reads in the Einstein frame as, 

^ +^^^^ =0, (17) 

where the semicolon represents a covariant derivative which may be expanded in terms of the 
Christoffel symbol as T^-^ = T^,^ + T^afj.T" — T^j^^Ta- An observer unaware of dark energy 
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looking out at the universe would come to the conclusion that energy is not conserved. That is 
not the case, however; it is the combination of matter and dark energy that is conserved when 
there are interactions in cosmology's dark sector. Summing the terms and remembering that 



Tu is diagonal and the relevant ChristofFel symbols are F^qq = 0, and f = f^Qg 
gives the conservation equation for a perfect fluid 



n 



03 



Pi + 3Hpi (1 + Wi 



0, 



H 



(18) 



where pi oc {aA{(j))) 3(i+Wj)_ The fractional energy densities are defined as 

so that for each component < Qi < 1. Using the expressions pi = A^pi and a = A{(f))a and 
carefully considering the derivatives involved, the fluid equation becomes 

^ dA (0) 



Pi + 3i7 (1 + Wj) Pi 



(1 - 3wi) Pi. 



A{(p) dcp 

Out of the four equations, (jl3p -(|15p and (j20p . only three are independent 



(20) 



55|, 



56( 1 . That is to 



say, Eq. (f20l) may have been derived from (fT3]) - (fT5]) . without assuming (fTTl) . 
From Eqs. (llSp and (I14p . we find that 



a 

a 



AttG 



{pi+3pi) + (p^+sp^) 



i=l 



(21) 



A species that contributes to positive cosmic acceleration (a > 0) must have (p + 3p) < 0, i.e. 
the equation of state w must be < —1/3. This gives us limits on the range of values the equation 
of state which the scalar field may have over time if it is to contribute a gravitationally repulsive 
force, < —1/3. Indeed, the WMAP data combined with BAO and SNIa observations has 
put much stronger limits on the DE equation of state, at present (z ~ 0), it lies in the range 
-1.11 < WDE < -0.86 (95% CL) H. This range includes values of equation of state less than 
— 1, which is outside a theoretical limit set for a canonical scalar field and belongs to "phantom 



cosmology". However, it is possible to get at 
to couple with dark matter non-minimally 



551, 



ow redshifts wj, < — 1 by allowing the scalar field 



56l |. especially by allowing the coupling between 
(j) and dark matter to grow with time, or allowing A{(j))'^ to increase monotonically with (p. 

It is practical to express the above equations in terms of the e-folding time, defined by 
N= \n[a{t)]. This is a useful time parameter when one is considering the expansion history of 
the universe, including late time cosmology. We also utilize the identity d(j)/dN = -jj^ and 
denote the differentiation with respect to by a prime, d/dN = , and the time derivative 
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represented by a dot, d/dt = Using following definitions 

H d\ii [A{(l))] _ Pi _ _ 



one can write the equations (jl3p - (|15p and (j20p in the following form 

"^ni + n^ = l, (23a) 

i 

2e + 3 (1 + w^) 0^ + 3 ^(1+Wi) !^i = 0, (23b) 

i 

n'^ + 2en^ + 3n^ (l + w<^) + (^'Q^(l-3wi)f^i = 0, (23c) 

i 

^J^^ + 2e^a + 3 (1 + Wi)-0'QX1(1"^'^*)^* = O' (23d) 

i i i i 

where the sum over i represents the sum over all forms of matter or energy. It is important to 
note that in minimal coupling case, A{(f)) = 1, Q vanishes and the field equations have one less 
parameter making them easier to handle. Solutions of these equations can be used to reconstruct 
the behavior of (j) over time. The behavior of the scalar field is generally characterized by its 
potential, which is given by 

(0) = H\^) (^Sn^ - ^ <p''^ . (24) 

For simplicity, one can make a specific ansatz for the form of the potential in order to reduce 
the number of degrees of freedom and hence solve the set of field equations analytically. For 
many potentials discussed in the literature, the initial value of (f) and (p must be finely tuned to 
obtain the correct values of fi^ and w^ today. The tuning of the initial field expectation value 
is required in addition to tuning the potential parameters. 

Given the many alternative form of (quintessential) potentials it is useful to try and under- 
stand the properties of DE in a model-independent manner. In this paper, we will make an 
ansatz for the form of the scalar field so that the form of the potential, together with a number 
of cosmological parameters, can be reconstructed. This is a valid way to study the effects of the 
scalar field on the background since it is required that the potential must be fairly flat today, 
or, equivalently, satisfy that {(p/mpi) ^ \/^H. This is because the scalar field must be rolling 
slowly (with a small kinetic term) to force the equation of state close to —1 at z > 0, as inferred 
by cosmological observations (the WMAP data combined with BAO and SNIa observations). 

An unambiguous expression for the quintessence equation of state (EoS) is given by 

w,(a) = — ^ . (25) 
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This result is obtained by substituting the expression for V{(j)), Eq. (j24p . into the definition of 
Wff, in Eq. (j22p . In this note, rather than parametrizing the dark energy equation of state 
with two or more phenomenological parameters, we will choose to parametrize the evolution of 
the field (p (cf see Eq. (j30p ). This, in turn, will induce a dependence of w,^ on those parameters 
through Eq. (j25p . The same applies to the scalar field potential V{(f)) and the Hubble parameter 
which both depend on these parameters through Eqs. (j24p and (|42p . In our analysis, we only 
consider a canonical quintessence with 7 = 1, and use the relation (|25p which is an exact and 
definitive expression for the equation of state. 

Other parameters of interest are the effective equation of state and the deceleration parame- 
ter. The effective equation of state describes the average equation of state, taking into account 
all the energy components of the universe: 

Weff = ^, Ptot=P^ + J^Ke^-Q"^, p,,,^p^ + ^p,e4-«"^. (26) 

Ptot 

From Eq. (j23bp and the definitions in ()22p . we find that 
Another important quantity is the deceleration parameter 

This parameter is defined such that it is negative for an accelerating expansion (a > 0) and 
positive for a decelerating expansion (a < 0). 



Weff 



2e 



(27) 



III. MOTIVATION AND ANALYTIC SOLUTIONS 



In order to be able to analytically solve the generalized field equations obtained in Sec. II, 
there must be more assumptions made to reduce the number of degrees of freedom. In this 
paper, instead of making an ansatz for the scalar field potential or a paramterization of the dark 
energy equation of state, we wish to reconstruct relevant cosmological parameters from some 
phenomenologically well-motivated Ansatze for a scalar field quintessence. Our approach is new 
in the existing literature. 

In the literature, basically, there are two simple and convenient approximations to the form 
of a time-evolving scalar field quintessence, see, e.g. Refs. 75l477l|. Those approximations are 



K\(t)\ = a^hit + const, = a^f^ + const, (29) 
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These evolutions for (p generally represent the tracking limits of some more general solutions for a 
runway dilaton in many string-inspired scalar-tensor theories {2ol. [tsI.ItgI. Similar approximations 



have been adopted a 



Refs. 



60 



61 



77 



so in a few phenomenologically motivated models for quintessence, see, e.g. 



81 



S3|. 



In this paper, we will assume that the time-evolution of quintessence scalar field (j) is well 
described by a combination of the above two specific solutions. Using a general feature of the 
scale factor, that a{t) oc t^, with different values of p at different epochs, such as p = 2/3 in a 
matter-dominated epoch and p > 1 in a dark-energy-dominated epoch, the time-evolution of (j) 
may be written as 

K{(j)Q- (t){a)) = a In a + jSa'^'^, (30) 

where a, f3 and C, are some new arbitrary constants, and a = a{t) is the scale factor of a 
four-dimensional FRW universe. This parameterization of (p, which might bear some generic 
features of a fundamental scalar dilaton field or metric moduli in some string-inspired models, 
can also be motivated from other two aspects. First, it really gives a useful information as 
regard the dark energy equation of state once the parameters like a and or their combination, 
are known (even approximately) using observational data. Second, it certainly helps to explain 
the cosmic coincidence problem, since the value of an arbitrary coefficient Vq that arises with 
a prior choice of the scalar field potentials, such as V{4>) = Vq e~^^'^^"^^\ ^{(p) = and 
V{<j)) = Vo (l — co</>^) e~^^'^^^'^^\ won't be very important in our approach. Indeed, the param- 
eterization (130p for the evolution of (p is neither more arbitrary nor more restrictive than other 



pa,rametrizations and approaches to quintessential dark energy in the literature. See Refs. 



8C- 



85( 1 for some other plausible ways of reconstructing dark energy or quintessential potentials. 
Our ansatz for (j), i.e. (|30p . has nonetheless some similarities with respect to numerous dark 



energy potentials proposed previously, which may be reconstructed by using Eq. ()24p . or alter- 
natively, V{(j)) = [(3 + e)(l - ^m) + ^^'m] H'^{4>)- Especially, with C in ([30]) , we find that 
the leading term of the reconstructed potential is simple exponential in 0, 

V {(!)) = e'^^'^'^P^ + ■■■ , (31) 

where the dots represent some other terms which could arise, for instance, due to the effects of 
matter-scalar interactions. Similarly, with a ~ 0, the scalar potential takes the form 

V{^)^V^e-^^" {l-c^ip^)+-- - , (32) 

where ip = \/S.4>\/mp = {(p^ — (p)/mp and Cq = A polynomial potential multiplied with 

exponential pre factor as above may be motivated by string theory and standard Kaluza-Klein 
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gravity. An effective dark energy potential as above was proposed in [86|]. In our approach, the 
total or effective quintessence potential is roughly given by a Unear combination of the above 
two potentials. Interestingly, in our discussions below, the actual form of the potential is not 
very important, but only the values of the slope parameters a and ^. 



591]. The model 



In fact, using the freedom to rescale N (= In a) or shift <f>, we can set (3 = 1 
then contains two free parameters, a and so there is more freedom to tune them according 
to the observational constraints. In our model, there is one more degree of freedom that must 
be constrained. This extra parameter, which is the coupling constant between the dark matter 
and the scalar field, is a mixed blessing - it makes the system easier to tune but it also makes 
the solutions more complicated. 

In order to completely solve Eqs. ()23ap - ()23dp . one must specify some initial conditions. These 
conditions can be defined on the current composition of the universe, as measured by WMAP 
plus other observations Q], and are — 0.05, 0,cdm — 0.22, ~ 10~^ and ~ 0.73. The 
baryon component (subscript b) is quite small compared to the cold dark matter component 
(subscript CDM), so even though there are tight limits on any baryon-dark-energy coupling, 
it is sufficient to assume a general dark energy coupling to fim, (f^m = + ^CDM ~ 0.27). 
Henceforth the only components of the universe that are not considered negligible are dark 
matter and the dark energy component, up to z ~ 0(10^), when radiation starts playing a 
major role. Also, it is assumed that the dark matter equation of state is the same as that of 
ordinary matter (or dust), w^ ~ 0. Using these assumptions, and for simplicity setting k = 1 
from now on, the system of equations ()23ati23d"|) is simplified to 

Q.m + ^4, = 1 (33a) 

2e + 3 (1 + w^) 0^ + 3f]„ = (33b) 

J^'^ + 2eO^ + 3f]0 (l + w</,) + 0'QO^ = O (33c) 

n'^ + 2en^ + ?,n^-(l)' Q^rn = Q- (33d) 

Equations ([33al) - p3dl) can now be solved to find explicit expressions for w^, e and Vtm- Solving 
the equations directly gives the solution for 0.^ and e in terms of a, Q and Q (a and Q are 
parameters from the ansatz for and Q is the coupling term). 
The explicit solutions for Q.m and e are given by 

""- C„-3JXW^^ ' -i,-.-2n„-i(„ + 2C.^C.)\ (34) 

where 

X = X{N = Ina) = exp [-(3 - aQ - a^)N + (2a + Q)e^^^ + Ce^^^l . (35) 



13 



The scale factor a{t) is normalized such that a{z = 0) = = 1, so in the past a < 1 and hence 
= In a < 0. The relationships among the scale factor a, redshift z and e- folding time N are 

— = -^, N = \n{a/ao) = -ln{l + z). (36) 

So, at t = t^, z = N = 0. Without loss of generality, we set oq = 1, so = In a. To solve 
the integral (I34p . in the small (" limit, we may use the approximation e^^'^'* w 1 + ^Ina. The 
solution for is now given by 

3 + Ciexp[(3-a2_Qa)Ar]' ^'^^> 
where we have introduced a new variable d such that 

d = a + 2C. (38) 

The solution for automatically gives the solution for il^ which is simply i}^ = 1 — il^, 
in using Priedmann constraint. The integration constant Ci can be fixed by using the initial 
conditions of the model, assuming r2m(io) = ^mO' where to is the time now. 

Particularly, at low redshifts, the terms quadratic (and higher powers) in (^N contribute only 
subdominantly. [Even at high redshifts, the terms like (2a + Q)e'^^'^ and ^e^''^ in (I35p are only 
sub-leading to the first term in X{N), i.e. — (3 — — aQ)N, since C > (by assumption) and 

= Ina < in the past.] To quantify this, one can write 

a2i„ „2 

gClna = i + <^lna+ ^ ^ +••• = 1 + Clna + C^|lna| +••• . (39) 

For instance, in between the redshifts z = 2 and z = 0, iV = Ina runs from —1.09 to 0, and the 
solution ()37p remains valid for as large as ^ ~ 1/2, but in the discussions below we will always 
assume rather implicitly that ( <^ 1/2. 

The Hubble parameter may be evaluated by solving the differential equation associated with 
the relation He = H' . This again introduces an integration constant which can be fixed by the 
normalization H{tQ) = Hq. From the solutions for and e, the other parameters discussed in 
Sec. II can also be derived. The expression for w^ is given by rearanging Eq. (j33bp . while V{4>) 
is given by Eq. ()24p using the solutions for flrn. and H. 

From the second equation in (|34p . we can see that (retaining the 4d gravitational constant 

^) 

^<p" = l + q-ln^. (40) 
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Requiring that q < (as implied by the type la supernovae) and imposing a generous lower 
bound on the value of il^, which is > 0.24, one obtains the safe upper bound 

k\cP'^\ < 0.8. 

However, to make the present model compatible with various other data sets, including WMAP 
observations, one may be required to satisfy K\(p'^\ < 0.4, which encompasses the cosmological 
constant limit within la error. 



IV. CONSTRAINTS ON THE MODEL 



A. Constraints from supernova 



In this subsection, we use several sets of data from recent cosmological observations and put 
constraints on our model, limiting the values of a, (" and Q. To do this, code for curve fitting 
given by Nesseris and Perivolaropoulos 72] has been utilized along with our solution for w^ and 



the standard form for the Hubble parameter, in terms of the two parameters 

pi = d^, p2 = -Qa, (41) 

where, as before, a = a + 2C,. We use the Hubble parameter to fit our model to data on the 
expansion rate of the universe. The Hubble parameter can be found from the solution for e and 
and the solution to the differential equation, e = H /H^, which is the definition of e. After 
some simplification, we find that the Hubble parameter is given by 

H = H^ ifloma-^-P^ + (1 - Oom) a-P^f'^ . (42) 

This result is obtained by integrating out the expression e = H' /H = (Ini?)' (see, Eq. (j34p ) and 
substituting the expression of Q.m from Eq. (j37p . The frame where matter density decreases as 
Pm oc l/(A((/>)a(t))'^ (where the scale factor is modified by the coupling A{(f))) is known as the 
Einstein frame. In our expression for the Hubble parameter (j42p the matter part is Oom^J^^"^^) 
where the density evolution is modified by the parameter p2 which is associated with the coupling 
Q. Thus this standard form of the Hubble parameter is in the Einstein frame. 

As in some standard approaches, let us first drop the scalar field - dark matter coupling. The 
model then reduces to one-parameter parameterization of the Hubble expansion rate. Remember 
that the model still deviates from the ACDM model, since Pi 7^ even though p2 = 0, leading 
to a time- varying equation of state for dark energy (see cf Fig. [T]). Needless to say, the ACDM 
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cosmology corresponds to the choice a = = Q. This can easily be understood from the 
following observation. With K(f)' ~ a + 2^ = a, we get 



_ «^ - 3Q^ _ ia^-3)C,-3aQ{l + zf-'' -"Q 

3C, + 3a{a + Q)(l + z)3-"^-«Q " ^^""^ 

in the above equation or in (j37p is fixed such that fiml^o) = ^mo at z = 0. Note that the 
dark energy EoS w^(2;) is time-varying as long as a 7^ 0. 

One notes that a prior choice or an arbitrary parameterization of w{z) could easily 
lead to an erroneous reconstruction of the dark energy equation of state. To see this, 
one considers the parameterization H'^{z) = //q [r2om(l + -2)^ + ^^de] , where J^de = (1 — 
riom)exp|3 J For example, with the ansatz w{z) = wq + j^wi, one finds that 

^DE = (1 — ^Om) exp I + (1 + Wq + W-,^) ln(l + z) — wi + constj. Now, with wi = 0, one 
has H'^iz) = Hi [J7om.(l + zf + (1 - r^om.) x const x (1 + 2;)3(i+"o)j ^ which has a form similar 
to Eq. (ji2]) . especially, with p2 = 0, i.e. (5 = 0. This could give a wrong impression that in our 
model the choice (5 = gives a constant equation of state for dark energy. The dark energy 
equation of state is not constant as long as a 7^ 0, or precisely, when a, C 7^ 0. 

In Table 1 we present the best-fit values using only one parameter (i.e. 7^ and p2 = 0) 
first using only the Gold sample of 157 type la supernova data, the Supernova Legacy Survey 
(SNLS) data alone and then for combined data sets. The combined data includes the cosmic 
microwave background (WMAP) shift, baryon acoustic oscillations (BAO), suvernovae type la 
Gold sample (SNIa) and legacy survey (SNLS). The errors of these fits are shown in Fig. [TJ 
For the combined data sets is minimum when f^om ^ 0.27 where a + 2<^ ~ 0.3 it 0.3. For 
^Om > 0.29 in the fit to the Gold sample alone, the EoS drops below -1, indicating phantom 
quintessence and an imaginary a. 

It is indeed the SNLS data that lowers w^{z) towards the value —1 at present, i.e. at z = 0, 
which is clearly seen from the best-fit values in Table 1. The SNLS data naively suggests a 
small cross over range between the cosmological constant w^ = — 1 and the phantom dark 
energy w^ < — 1 (or, equivalently, pi < 0), but the error bars are too large. For the combined 
data sets (from WMAP, BAO, SNIa and SNLS) the best-fit value of w,^ falls in the range 
-0.94 > w,^ > -1. 

From the relationships above, and also noticing that for the combined data sets the scalar 
field equation of state is closer to —1 over a longer time in Fig. [U we find that the solution 
using all the data sets is compatible with a cosmological constant. This analysis has provided a 
limit to the relationship between a and C, but there is still a degree of freedom for choosing the 
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Table 1: The best fit of a to expansion history data, Q = 0: 



Om 


V\ 


|a| 


w^{z = 0) (mean) 


mm 


SNIa Gold data sets 


0.26 


0.271 ±0.286 


0.521 (±0.275) 


-0.878 


177.98 


0.27 


0.194 ±0.298 


0.441 (±0.339) 


-0.911 


177.76 


0.28 


0.112 ±0.312 


0.335 (±0.466) 


-0.948 


177.54 


SNLS data sets 


0.26 


-0.006 ±0.274 




-1.00 


104.15 


0.27 


-0.081 ±0.285 




-1.04 


104.13 


0.28 


-0.159 ±0.297 




-1.07 


104.12 


WMAP+BAO+SNIa+ 


SNLS data sets 






0.26 


0.123 ±0.178 


0.350 ±0.254 


-0.945 


283.54 


0.27 


0.098 ±0.1842 


0.313 ±0.294 


-0.955 


283.24 


0.28 


0.072 ±0.190 


0.268 ±0.354 


-0.967 


283.61 



-1.5- 



0.5 



1.0 



1.5 





FIG. 1: Scalar field equation of state as a function of redshift (z) with Icr and 2ct errors (light and 
dark gray shades) and with EoS solutions given by C = 0.2, 0.15, 0.1, (top to bottom; green, pink, red, 
blue online), 17o„. = 0.27. Left plot: SNla+WMAP+BAO+SNLS. Right plot: SNla only 



value of a or ^. This degree of freedom proves hard to constrain since for other observational 
tests, it is only required that a and C, be of comparable magnitude, or that C = 0. The plots 
of W0 in Figure [1] show that the chi-squared is minimized when C, is small since cosmological 
reconstructions for large ^ diverge from the best-fit line. Figj2] displays the same constraints on 
the model, this time showing the effect of positive or negative C,. There is more divergence from 
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the best-fit line for negative values of C, suggesting that positive is a better fit. 



0.5 



1.0 



1.5 




-1.5 




FIG. 2: As in Fig. [T]but with = —0.2 (lower line, eyan) and ( = +0.2 (upper line, blue). Left plot: 
SNla+WMAP+BAO+SNLS. Right plot: SNla only 
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-1.0 



-1.5 
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1.0 y 


/ 1.5 


= 













FIG. 3: Best-fit plot of w^ with Q (i.e. p2 7^ 0). The continuous line is the best fit and the dashed 
line is for ~ 0.2. The inner and outer shaded regions represent Icr and 2a error bars. Left plot: 
= 0.26. Right plot: l^o™ = 0.28. 



Our proposed model has the freedom of allowing for a non-zero coupling of dark matter to 
dark energy. Including this in the fit to the data sets gives the results of Table 2. Note that 
for Qq > 0.27 we get negative values of pi. From the definition of pi, this means an imaginary 
a and Q, and a negative equation of state which is phantom quintessence. The best-fit to the 
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Table 2: The best fit values for both a and Q (WMAP+BAO+SNIa+SNLS): 



Om 


Pi 


P2 


|a| 


Q 


w^{z = 0) 


mill 


0.26 


0.055 ±0.203 


0.015 ±0.021 


0.235 ±0.432 


-0.065 ±0.149 


-0.98 


283.02 


0.27 


0.001 ±0.210 


0.021 ±0.021 


0.027 ±3.83 


-0.772 ± 1.08 


-1.00 


282.22 


0.28 


-0.055 ±0.217 


0.027 ±0.021 






-1.03 


281.98 


0.29 


-0.111 ±0.225 


0.032 ±0.021 






-1.05 


282.24 


0.30 


-0.168 ±0.233 


0.037 ±0.021 






-1.08 


283.00 



data is for il,„o — 0.275 which is apparently the best fit value for ACDM model obtained from 
the combined WMAP5+BA0+SNIa datasets The best fit value for a is close to zero, the 
cosmological constant limit. Q is constrained to — 0.772 ± 1.08 but this error is too large for this 
to be a conclusive result. It is clear from the shaded error region in Fig. [3] that adding the extra 
parameter, Q, greatly increases the uncertainty in the fit. 



B. Constraints from structure formation 



c energy on the expansion 



63l |. Other than this there 



Cosmological models are probed by observing the effects of darl 
of the universe through measuring the distances to far off galaxies 
are not many ways to check models that account for the late time acceleration of the universe. 
However, if dark energy does interact with dark matter even weakly, it may have an observable 
effect on the early stages of structure formation. 

Since dark matter naturally guides the way for the formation of observed structure, any 
interaction between dark energy and dark matter would have an effect on the manner that 
visible structure was formed. In particular, we investigate how the rate of structure growth is 
affected by a non-zero coupling of dark matter to a fundamental scalar field. 

Matter fluctuations evolve according to the standard linear differential equation 

5 + 2H5 = ATrGpm5, (44) 
where 5 = 5pm/Pm is the linear matter density contrast. This linear growth equation comes 



M\. It 



from the perturbed equations of motion of Einstein's general relativity, see, e.g. Ref. 
may be numerically solved to reconstruct the growth of these matter fluctuations, but as yet 
there is no physical theory that relates the matter density contrast to the matter density. In 



65( it was proposed that the growth rate of matter perturbations, defined as 
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can be characterized by the following simple expression 



f=[nrn{N)]\ 



(46) 



This ansatz works well at low redshifts {z < 2) and for c oup led dark energy models with a 



65( 1 ■ and for modified DGP gravity, 



small coupling parameter. For the ACDM model, r] ~ 0.56 
r] ~ 0.68 661]. However, this ansatz does have the drawback that / < 1, since < dm ^ 1- This 
is fine for the standard ACDM model which assumes / = 1 during the matter-dominated era 
for high redshift, but does not hold for all models of dark energy. 

For the present model of quintessential dark energy, as discussed in 
linear growth equation (j44|) is modified to be 



671 ] in more detail, the 



6'^ + {2 + 6 + Qc (P') 'Jc = ^ (1 + ScQc + ^(1 + ^QMSb^b 
6'^ + {2 + 8 + Qb(t>') 5'b = l{l + 2Ql) dbQb + ^(1 + 2QbQc)SA 



(47a) 
(47b) 



for baryonic and CDM components respectively. The couplings constant Qb and are usually 
coupled, so in g eneral it is not possible to express these equations as a single differential equation. 
Amendola 



671 ] made a naive estimation that in the limit \Qb\ <^ \Qc\ and ^Ib ^ one can 



write 



6" + {2 + e + Q (t)')5' = + 2Q^)6, 



(48) 



where Q = and 5 ^ 6c (to leading order). In a sense, the baryonic component is assumed to 
be negligible. But one should also note that 6b and its derivatives are non-negligible, otherwise 
Eq. ()47bp would be inconsistent. Equation (|48p . which reduces to the standard result when 
Q = 0, is a close approximation rather than being an exact result. In the following discussion, 
we will assume that 6c^c ^ 6b^b- 

In space-time backgrounds dominated by baryonic matter, the effect of dark matter may be 
neglected {^c = 0). The effective Newton's constant G, and two dimensionless post Newtonian 
parameters 7 and /3 are related to the coupling constants 



via 



Q{(f)) = mp 



G = G [A\cj)){l + Q"" 



din A{(l)) 



X{(j)) = mp 



dQ{^) 



'<^o' ^ - 1 + Q2' ^ - ' 2(1 + Q2)2- 
For an exponential coupling ^(</>) oc e^'''^'^, with Qb behaving (almost) as a constant, the local 



1 



2Q2 



/3 = 1 + 



(49) 



(50) 



gravity constraint ]1 — 7] < 10 



691] implies that = {1 - j)/{l + 7) < lO"'' and hence 
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z 


0.15 


0.35 


0.55 


0.77 


1.4 


3 


fobs 


0.51 ±0.11 


0.70 ±0.18 


0.75 ±0.18 


0.91 ±0.36 


0.90 ±0.24 


1.46 ±0.29 



IQbl < 0.01. This constraint is still weaker than the one arising from weak equivalence principle 



violation 



71[. With Qb — const, one has |1 — /3| ~ to a large accuracy. The effective 
Newton's constant is now modified but it must be within the limit for the time variation of 
Newton's constant | '^'^'^^ | < lO^-'^'^yrs"-'^, see e.g. Refs 



HQ , which is generally the case when 
\Qb\ ^ 0.01 is satisfied. In most of the discussion below we will assume that \Qc\ ^ \Qb\ ~ 
and Qc = Q- 

In the nonminimally coupled theory, the linear growth rate (144p is modified from the uncou- 
pled case; the quantity 6 = 5pm/Pm now depends on the values of Q and (f)'. The result of this is 
a modification of the expression for the growth rate, Eq. (|46p . where it is renormalized by either 
the parameter 77 — )• ry and/or a coefficient /q: 

/ = /o [nmiN)]" . (51) 

Following the analysis in Ref. [tJ], if the growth rate is modified by a coefficient, it may be 
dependent on the coupling by 



(52) 



which reduces to (j46p for minimal coupling (Q = 0) . This expression of / is merely a phenomeno- 
logical fit where rj and c may be determined by fitting the standard approximate solution (j52p 
to a numerical solution to Eq. (j47ap or Eq. (08]). The best-fit parameters found from this nu- 
merical analysis in Ref. [tJ] are 7] ~ 0.56 and c ~ 2.1. The physical theory pertaining to the 
dependence of and Q on the perturbation growth rate is as yet unknown so this kind of 
generalized fit is just an approximation. Determining the growth rate over time and comparing 
to observational results gives limits to the model that are independent of limits from fitting the 
model to supernova and WMAP data sets. This grants an independent check for the consistency 
of the model. 

To preform this check of our model, we use some known observed values for the growth rate 
at different low redshifts (listed below in table), from which we can determine the values of /o 
and r] that best match the data. There are very few measurements of the growth rate at low 
redshifts and often / can be evaluated by assuming a ACDM model, for which r] ~ 0.56. 

It is clear from Eqs. (j37p that Qm is in fact a function of a and Q and thus our parameteri- 
zation of /, i.e. ()52p . is implicitly a function of a. The fits are illustrated by Fig. IU where the 
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best-fit parameters are 

\Q\ = 0.199, a = -0.010, r/ = 0.538. (53) 

In the above we have used the normahzation /o = 1 + 2.1Q^ used in {tJ]. Within the la errors, 
the normalization constant /o = 1.16 it 0.34 is compatible with a minimally coupled quintessence 
scenario. The least-squares curve fitting of all three of a, Q and t] to the data often leads to 
huge ranges of possible values of these parameters. For brevity, we may restrict ourselves to 
fitting only two parameters such as Q and c. The best-fit value of the normalization constant 
c is quite sensitive to the choice of or vice versa. For example, assuming QmO = 0.27, a = 
and T] = 0.56 and fitting for Q and c gives 



Q = 0.42 ±0.10, c= 1.17 ±0.02. 



Similarly, holding r] = 0.68 and c = 2.1 give the values 



IQI = 0.41 ±0.14, d = 0.57 ±0.38. 



(54) 



(55) 



In the above analysis we have used an explicit analytical expression for which is nothing 
but the equation (I37p . If it is instead the parameter ij that is renormalized, say rj ~ 0.68 66 1 
as in modified gravity proposed by Dvali, Gabadadze and Porrati, the fit to the data in Fig. |4] 
seems to be closer. 




FIG. 4: Best fit curve to the observed values of the growth rate, with fiom = 0.27. The red solid line 
eorresponds to a best fit of Q, a and 77, the blue dashed line corresponds to 77 = 0.56 with best fit for Q 
and a, and the blue dotted line corresponds to 77 = 0.68 with best fit for Q and a. 



C. Constraints from CMB 

Another test of our model to the observational constraints is to demand that the equation of 
state presently lies in the range, —1.11 < w^ < —0.86, as implied by the WMAP data combined 
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with BAO and type la SN Q]- This is illustrated in the contour plot, Fig. [5l where the outer 
contour corresponds to w,^ = —0.91 and the innermost contour corresponds to w^ = —0.99. For 
w^ to lie within the WMAP5 limit —1.11 < < —0.86, the constraint is 



a = 0(±0.5) 



(56) 



where the value of the coupling has little impact. This relationship holds for small \C,\ < 0.1. 
It is reasonable to assume that the value of the scalar field equation of state has not changed 
much in the recent past. Using this tighter restriction, by inspecting Fig. [5l we get a tighter 
relationship between a and as z — t- 3, error in relationship (|56|) goes to ±0.2. 




0.1 0.2 0.3 0.4 



OC -1.0 




FIG. 5: Left plot Contour plot of in the range (—0.91, —0.99) (top to bottom) at z = 0. Right plot 
The dark energy EoS for the two solutions of C = ±0.05 (solid lines) and C = ±0.25 (dashed lines). 



There is still a degree of freedom in the choice of the value of a or C, so we look to a possible 
link between inflation and late time acceleration to set it. Quintessence is designed such that 
the scalar field that drives the late time acceleration may be the same, or an evolved form of 
the infiaton field that drove infiation. If we assume that it is in fact the same field that drives 
both of these accelerating epochs, we may impose extra conditions on the model by requiring 
the quintessence field to satisfy the constraints of infiation. One such constraining parameter is 
the spectral index, which describes the slope of the angular power spectrum of the CMB. The 
WMAP data inferred a red-tilted spectrum, < 1, which is consistent with most inflationary 
models. The spectral index can be approximated in the small, positive C, limit, as 4^ 

~ 1 — a^. (57) 

For the WMAP±BAO±SN mean value of = 0.960 l|5;2i3 and |a| = 0.20 ± 0.04. By using 
the relationship between a and C, determined by the combined data sets in Table 2 (|q ± 2C,\ ~ 
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0.3), we obtain ( = ±0.05 ± 0.15 or C = ±0.25 ± 0.15. The evolution of the EoS for these 
solutions of C are plotted in Figure [5j The EoS for ( = ±0.25 deviates far from -1 for a positive 
value of a. Both signs of a are possible so the lower value C = 0.05 may be more physical. 

V. DYNAMICAL BEHAVIOUR 

In this section we examine how the scalar field varies over time and how this depends on the 
parameters a, ( and Q. As an example we choose the value a = 0.2 and using a + 2C ~ 0.3 from 
Table 1, C = 0.05. From structure formation we estimate \Q\ = 0.2. In our analysis perhaps 
the most relevant parameter is the field velocity of <j), i.e. (j)' = d(j)/d\na, which is a measurable 
quantity, at least in principle, either by accurately measuring the equation of state w^{z), or by 
placing a constraint on the dark matter-scalar coupling Q, or both. 




FIG. 6: (Left) Form of the scalar potential implied by ([M)) . with Q ~ 0.2,0 and —0.2 (top to bottom) 
(blue, green and red). (Right plot) Shape of a reconstructed scalar field potential. 

One can think of the scalar field "rolling down" the potential, like a skier on a mountain. It 
can be seen from Fig. [6] that the slope of the potential now is not quite flat; rather, the scalar 
field is rolling with a changing velocity. This is apparent when taking the time derivative of 
the scalar field, (j)' = —a — 2C,e^^^ . At = 0, \(f)'\ = a + 2C,. Indeed, we require that the 
scalar field is rolling slowly so that the kinetic energy is small compared to the potential energy. 
The scalar field rolls at a slowly changing velocity for a small value of C, and the acceleration is 
(j)"{N = 0) = —4^^, so we expect that the value of to be much less than 1. Furthermore, cp" 
is negative so the scalar field is slowing down, it is in a 'freezing' phase. The model approaches 
the ACDM model as a, C 0. 

The shape of the reconstructed potential shown in Fig. [6] is essentially an exponential. In 
particular, in the limit —t- 0, it is a simple exponential, V{(l)) oc e°"^, and for the non-zero ( case 
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this is modified by extra terms proportional to e'^^ (0 — 4>of'. The reconstructed potential shown 
on the left panel of Fig. [6] displays a range of possible values of Q; for larger Q the potential is 
steeper, and for smaller Q it is flatter and even changes shape for a negative Q. 




FIG. 7: Scalar field equation of state as a function of redshift. We take the values as a = 0.2, C, = 0.05 
and ilom = 0.27. 

The equation of state for our choice of the parameters is currently within the limits from 
WMAP, —1.11 < W0 < —0.86, and does not change much up to redshift z ~ 1. The effect of a 
nonzero coupling Q on the equation of state is shown in Fig. [71 At low redshift, a positive value 
of Q drives w^ closer to —1, whereas a negative Q gives a high value of w^ at high redshift. 



VI. EFFECT ON THE BACKGROUND 



An intrinsic property of dark energy is that is does not interact with light. The only way 
that it may be observed is through its effect on the evolution of the background and its possible 
interaction with dark matter. In this section we look into how this model for quintessential 
dark energy affects the background by quantitatively looking at its effects on the rate of cosmic 
expansion and the fractional densities of dark matter and dark energy. A useful parameter used 
when considering cosmic acceleration is the "deceleration parameter", q. In an accelerating 
epoch, g < 0. In Fig. [HI we have shown a typical variation in the value of q with redshift. The 
important feature of this reconstruction is that q drops from positive to negative at z ~ 0.7, 
when cosmic acceleration began 



90t |. A nonzero coupling has an effect on the redshift when 
cosmic acceleration started, Zacc] the positive coupling [Q > 0) implies an earlier start, and the 
negative coupling [Q < 0) implies a later start. 
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FIG. 8: (Left) Deceleration parameter for varying a {( = 0.05, Q = 0.2). (Right) History of cosmic 
acceleration for Q ~ 0.2,0 and —0.2 (bottom to top) {( — 0.05, a = 0.2). 

An alternative way of looking at the start of cosmic acceleration is to see when the effective 
equation of state drops below —1/3 in Fig. [9l At high redshift, the effective equation of state 
goes towards 0, which is expected in the matter-dominated epoch (w^ — 0). For all values of 
a, the effective equation of state goes to ~ during the matter dominated epoch, but for large 
values of a, w^s > 0. This is because the relative energy density of dark matter and dark energy 
now is fixed and the matter equation of state is constant, but the scalar equation of state is 
closer to for a larger a, as shown in Fig. [5l At least, at low redshifts {z < 1), the effective 
equation of state does not vary much for different values of the coupling. 



W0 Weff 




FIG. 9: (Left plot) The scalar field equation of state and (right plot) the effective equation of state with 
C = 0.05, a = 0.2 and Q = -0.2, and +0.2 (from top to bottom). 

The relative energy density of dark matter to dark energy is affected by changes in the a 
parameter, as seen in Fig. [TOj For smaller values of a, the matter density tends towards a 
lower value at high redshift. That is, due to the Priedmann constraint (|33ap . there is a higher 
proportion of dark energy during the matter dominated epoch for a small a component. Here, 
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it is clear that there has been an assumption made for the matter density today. The right 
hand side plot of Fig. [10] illustrates the effect of a non-minimal coupling Q on the evolution 
of matter density. The effect is not great within this small range, but there is a difference. A 
smaller value of the coupling parameter, Q, leads to an increase in the proportion of matter 
in the matter-dominated epoch. This study is valid for a range of low redshift, but will break 
down at high redshift when we enter the radiation-dominated epoch because we have made the 
assumption that the radiation component is negligible. 
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FIG. 10: (Left plot) Matter density for varying a (C = 0.05, Q ^ 0.2) with rJom = 0.27. (Right plot) 
Matter density for Q = —0.2, and 0.2 (top to bottom; online: red, green and blue) (C — 0.05, a = 0.2). 



One requirement of the potential is that it should dominate the kinetic term so that the 
scalar field equation of state, w,^, may be driven towards —1. It can be seen that this is satisfied 
by examining Figure [9l Recently, the authors of [62| have found two new diagnostic tools for 
distinguishing quintessential dark energy from Einstein's cosmological constant. Their method 
may be extended to a non-minimally coupled theory for which \Q\ > 0. Here we simply note 
that for small values of a and (, satisfying the constraints given in table 2, the constraints of 
the present model encompass the cosmological constant limit at the la level at low redshifts. 



VII. CONCLUSION 



In this paper, a scalar-tensor theory of quintessential dark energy has been reconstructed. 
This reconstruction rests on the general action ([6|) , in which the tensor theory of general relativity 
is modified by including a fundamental scalar field. The scalar field affects the gravitational part 
of the action by introducing a gravitationally repulsive term which drives the cosmic acceleration, 
and is coupled to the matter part of the action, giving rise to a nontrivial dark matter - dark 
energy interaction. As for all quintessence models, this action reduces to the one of general 
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relativity, when the scalar field is time-independent and the coupling of dark matter to a scalar 
field (or dark energy) is minimal. 

We made an ansatz for the form of the scalar field itself, in order to be able to reconstruct 
relevant cosmological parameters, including ^}m,^4>,^ and w^. The ansatz \(j)Q — (pl/nip = 
a In a + exp (2(^ In a) + const for the scalar field was motivated by generic solutions of the dilaton 
field in some effective string theory models in four dimensions. The two new parameters, a and C 
as well as the matter scalar coupling parameter, Q, were constrained using several cosmological 
methods. Primarily, the relationship between a and C was constrained by curve-fitting to various 
cosmological data sets. For example, with il-om = 0.27 and minimal coupling, using just the 
SNla data, minimizing gave a + 2^ = ±0.44 (±0.34) and using all the data sets available 
gave a±2C = ±0.3 (±0.3). 

We also considered the effect of dark energy on the growth of large scale structure. The effect 
of a nonminimal coupling on the standard expression for the rate of growth of matter perturba- 
tions, / = [ilm(-/V)]^, was considered. A modified ansatz for the dependence of matter density 
and dark matter-dark energy coupling on the matter fluctuation growth rate was introduced 
and then fitted to some observed values of the growth rate at low redshift. This analysis gave 
a best-fit minimum value for the coupling Q = 0.2 ± 0.2 (with the input ilom — 0.27), which 
is compatible with the minimally coupled scalar field case. These methods of constraining the 
variables of the model have succeeded in constraining a ± 2^, but failed in constraining the 
separate values of a and C. This is a direct result of taking the small C limit when analyzing 
the solution for flrn.. Nevertheless, by assuming the same field for the inflaton as quintessence 
the individual values of a and ^ were constrained using the spectral index of the CMB. We find 
a = 0.20 ± 0.04 and C = 0.04 ± 0.15. 

Our analysis showed that the present model of coupled quintessence is compatible with a 
cosmological constant (where a, C = and Q = 0), within la error bars. This means that dark 
energy may still simply be a cosmological constant, but even a small deviation away from this 
limit results in a very different source of cosmic acceleration, with time varying behavior. The 
number of parameters of this model may make it hard to rule out, as they may be tuned to 
satisfy many more constraints. In this work, the matter scalar coupling parameter was chosen 
to be a constant. This was not motivated physically but was chosen for the solvability of the 
system. Indeed, there is no physical reason why the dark energy - (dark) matter interaction 
should not vary in time, it is in all likelihood a function of (f). There needs to be further study to 
explore this (j) dependence of Q. Also, the links between quintessential dark energy and inflation 
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deserve further investigation. 
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